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Abstract 

We have measured production cross sections and longitudinal momentum distri- 
butions of fragments from neutron-deficient 36 Ar at 1.05 ^4GeV. The production 
cross-sections show excellent agreement with the predictions of the semiempirical 
formula EPAX. We have compared these results, involving extremely neutron de- 
ficient nuclei, with model calculations to extract information about the response of 
these models close to the driplines. The longitudinal momentum distributions have 
also been extracted and are compared with the Goldhaber and Morrissey system- 
atics. 
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1 Introduction 



Projectile fragmentation at high energies (~ 1 A GeV) has proven to be a 
powerful tool to produce clean secondary beams of exotic nuclei that can 
be used for a wide variety of nuclear-structure, nuclear-reaction, or nuclear- 
astrophysics studies. To make full use of this potential of projectile fragmenta- 
tion, the reaction mechanism has to be understood in detail in order to predict 
precisely the properties of the secondary beams. Pioneering experiments were 
performed at the LBL Bevalac accelerator [1,2]. Later, more detailed studies 
have been undertaken at the SIS/FRS exotic beam facility at GSI, Darmstadt 
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[3,4,5,6]. Comprehensive datasets have been obtained with respect to forma- 
tion cross sections of exotic fragments as well as their kinematic properties 
(mainly centroids and widths of longitudinal momentum distributions). Both 
types of observables could be fitted by empirical parametrizations: the former 
by the widely used EPAX parametrization [7], the latter by the Morrissey 
systematics [8]. These parametrizations have turned out to be useful within 
simulation programs of projectile fragment separators (e.g. MOCADI [9] or 
LISE [10]). More insight into the underlying physical processes was provided 
by comparison to physical models like abrasion-ablation (ABRABLA [11,12]) 
or intranuclear-cascade models (ISABEL [13]). 

Since most projectile-fragmentation experiments aim at producing exotic nu- 
clei near the limits of the presently known regions, and since the fragment 
neutron- or proton-excess is correlated with the neutron- or proton-excess 
of the projectile [7], fragmentation of stable projectiles outside the valley of 
/9-stability has received most interest. Quite a few studies of neutron-rich pro- 
jectiles like 48 Ca, 86 Kr, and 136 Xe have been reported in the literature [2,4,14]. 
On the proton-rich side, fragmentation cross sections for 58 Ni, 112 Sn, and 
124 Xe have been published [3,15,16]. It was found that the empirical EPAX 
parametrization could fit the measured cross sections with good accuracy down 
to sub-nanobarn levels [7]. 

The present paper aims at complementing the data base for the fragmentation 
of proton-rich projectiles by studying 36 Ar fragmentation. The isotope 36 Ar is 
important for producing light proton-rich nuclei at the proton drip line; their 
properties can be compared to those from 40 Ar fragmentation [1]. In addition 
to the formation cross sections (down to a level of /ibarn), widths of longi- 
tudinal momentum distributions could be measured; both will be compared 
in the following to the respective systematics as well as to physical model 
predictions. 



2 Experimental procedures 

The results reported in the present paper were obtained as a by-product of a 
secondary-reaction experiment aimed at measuring total interaction cross sec- 
tions of Ar and CI isotopes [17] as well as measuring longitudinal momentum 
distributions of knock-out products in coincidence with 7-rays [18]. This is the 
reason for certain experimental limitations of the present experiment that will 
be discussed below. The experiment was undertaken at the SIS/FRS facility 
at GSI in Darmstadt, Germany [19]. SIS delivered a 36 Ar beam with 1050 A 
MeV in spills of 8 s duration with a repetition rate of 1/16 s -1 . Typical beam 
intensities varied between 10 8 and 10 9 ions per spill, depending on the cross 
section to be measured. The primary-beam intensities were measured with 
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Fig. 1. A schematic view of the FRagment Separator (FRS) with the detection set-up. 
The isotope identification was possible in both sections of the spectrometer by mea- 
suring the time of flight (TOF) between scintillators (SC) and the energy deposition 
in ionization chambers. Several position-sensitive detectors (TPC) allowed tracking 
of projectiles and fragments and momentum measurements of the fragments. 



a secondary-electron transmission monitor (SEETRAM) [20] that was cali- 
brated with respect to a fast scintillator detector. Fragments were produced 
in a thick (1625 mg/cm 2 ) 9 Be target located at the entrance of the FRS. 

Fragments were identified with respect to their nuclear charge, Z, and their 
mass-over-charge ratio, A/Z, by a combination of energy-deposition, magnetic 
rigidity and time-of-flight (ToF) measurements. The set-up used is shown in 
Fig.l. While in principle both observables can be measured in both sections 
of the FRS, only the first half up to F2 was used to identify and count the 
fragments. To this end, the energy deposition in an ionization chamber at F2 
was measured together with the time of flight (ToFl) between scintillators 
"SCI" and "SC2". Positions and angles at F2 could be determined with the 
help of two position-sensitive detectors (TPC). The ToF was calibrated with 
the primary 36 Ar beam from SIS with several known energies which allowed 
to determine the length of the flight path and the offset. Together with mag- 
netic rigidity (Bp) determined from Hall-probe measurements and corrections 
for deviations from the optical axis, the ToF gave A/Z of the fragments. At 
the same time, the velocity-dependent energy deposition of the beam was 
recorded and calibrated in terms of nuclear charge, Z. Typical results of the 
isotope identification obtained are visualized in Fig.2. The non-observation of 
the unbound nuclei 19 Na and 16 F allowed to verify independently the isotope 
identification. 
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Fig. 2. An example of the fragment identification in the first half of the FRS. The 
isotope identification is confirmed by the non-observation of the unbound nuclei 19 Na 
and 16 F. 



3 Data analysis and results 



In the present experiment, the requirements of the main secondary-reaction 
study dictated an experimental setup that was not optimum for the present 
purpose. The main limitation was the width of the Fl start scintillator ("SCI" 
in Fig.l) of only 2 cm in x- (bending) direction. This leads to important cuts in 
the fragment longitudinal momentum distributions such that several Bp set- 
tings had to be analyzed to cover the entire phase space distribution of each 
fragment. Each Bp setting allowed to produce a two-dimensional spectrum as 
the one shown in Fig.2. For each isotope selected by a two-dimensional condi- 
tion in Fig.2, its x-distribution in the F2 momentum-dispersive focal plane was 
converted to a longitudinal momentum distribution and properly normalized 
with the help of the SEETRAM beam current detector. The reconstructed 
magnetic rigidity distribution of the isotope 14 is shown in Fig. 3. 

The reconstructed Gaussians for each isotope yield directly the average longi- 
tudinal momentum in the laboratory frame and the momentum width. From 
their areas, production cross sections can be determined once the necessary 
corrections have been applied. 
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Fig. 3. Magnetic-rigidity distribution of the fragment O reconstructed from several 
normalized Bp settings. The envelope of the distribution is fitted to a Gaussian. 

3.1 Production cross-sections 

An important ingredient in determining the production cross sections is the 
transmission calculation. In bending direction, a Bp scan allows to cover the 
entire longitudinal momentum distribution of each fragment or at least to 
cover enough momentum space that a Gaussian fit to the data is possible. 
The same is not possible, however, perpendicular to the bending plane. This 
means that the angular acceptance for each isotope has to be calculated nu- 
merically, based on its measured longitudinal momentum width and assuming 
that longitudinal and transverse momentum widths are identical. The nu- 
merical procedures have been published in Ref. [21]. It is well known that 
the momentum widths depend approximately on the square root of the mass 
difference between projectile and fragment (e.g. Ref. [8]). Consequently, the 
transmission of fragments close to the 36 Ar projectile (like the Ar or CI frag- 
ments in Table 1) are transmitted with about 100%, whereas the isotopes 
12 ' 13 N are transmitted only with about 33% efficiency 

Further corrections to the isotope yields are due to the dead time of the data 
acquisition (approximately between 50% and 70%) and secondary reactions 
of the produced residues in the target (of about 10%). Finally, the observed 
number of counts of each isotope is converted to a cross section using the num- 
ber of target atoms and the number of incident ions as determined from the 
SEETRAM beam intensity monitor. The resulting production cross sections 
for 38 isotopes are listed in Table 1 and visualized in Fig.4. 
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Fig. 4. Production cross sections of neutron- deficient ib Ar fragments measured at 
1050 A MeV. Solid lines link EPAX [7] predictions. Dotted and dashed lines depict 
model calculations with the ISABEL intranuclear- cascade [13] and the ABRABLA 
abrasion- ablation model [11], respectively. 
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Table 1 

Production cross sections of neutron- deficient fragments from 36 Ar measured at 1050 
A MeV 



Isotope 


O~prod (b) 


Isotope 


O~prod (b) 


12 N 


(2.97 ±0.09) • 10~ 4 


22 Al 


(2.30 ± 0.30) • 10~ 6 


13 N 


(3.72 ± 1.19) • 10~ 3 


23 Al 


(2.10 ±0.30) • 10~ 5 


13Q 


(4.20 ± 0.20) • 10^ 5 


24 Al 


(3.34 ±0.13) • 10~ 4 


14Q 


(4.12 ±0.13) • 10~ 4 


25 Al 


(3.48 ± 0.52) • 10~ 3 


15Q 


(9.82 ± 0.29) • 10~ 3 


24 Si 


(2.00 ± 1.00) • 10~ 6 


16Q 


(3.76 ± 0.64) • 10~ 2 


25 Si 


(3.40 ±0.10) • 10~ 5 


17p 


(2.99 ± 0.09) • 10~ 3 


26 Si 


(3.52 ±0.11) • 10~ 4 


18p 


(7.58 ± 1.35) • 10^ 3 


27 Si 


(6.13 ± 1.04) • 10~ 3 


17 Ne 


(2.30 ±0.10) • 10~ 5 


27p 


(7.40 ± 1.40) • 10~ 6 


18 Ne 


(5.67 ±0.17) • 10~ 4 


28p 


(2.33 ± 0.09) • 10~ 4 


19 Ne 


(3.12 ±0.09) • 10~ 3 


29p 


(2.78 ±0.17) • 10" 3 


20 Ne 


(1.80 ±0.40) • 10~ 2 


29 s 


(3.20 ±0.10) • 10" 5 


20 Na 


(5.24 ±0.15) • 10^ 4 


30 s 


(2.73 ±0.11) • 10~ 4 


21 Na 


(2.63 ± 0.38) • 10~ 3 


31 S 


(4.17 ±0.39)) • 10~ 3 


22 Na 


(1.47 ±0.46) • 10~ 2 


31 CI 


(1.20 ±0.30) • 10~ 5 


20 Mg 


(1.60 ±0.30) • 10~ 6 


32 CI 


(4.36 ± 1.32) • 10~ 4 


21 Mg 


(2.80 ± 0.50) • 10~ 5 


33 C1 


(6.59 ± 1.87) • 10~ 3 


22 Mg 


(5.19 ±0.21) • 10^ 4 


33 Ar 


(3.40 ±0.10) • 10~ 5 


23 Mg 


(5.13 ±0.15) • 10~ 3 


34 Ar 


(5.73 ±0.19) • 10~ 4 



The errors involved in the cross section determination are partly constant 
(2% each for the angular acceptance determination [21] and beam intensity 
determination) and partly vary appreciably from isotope to isotope. The main 
contribution comes from the errors of the Gaussian fit to the reconstructed 
momentum distributions. The total combined errors of the production cross 
sections vary between 3% and about 30% with the exception of 24 Si where the 
error is of about 50%. 
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Fig. 5. Widths of the longitudinal momentum distributions plotted as a function of 
mass loss, AA, from the projectile. The dashed curve shows the \J AA dependence 
of the Morrissey systematics [8]. The dotted curve follows Goldhaber's model [25]. 
The full curve has been obtained by fitting the constant <To in Goldhaber's model to 
the present data (see text). 



3.2 Momentum distributions 



As mentioned earlier, the momentum distributions for each isotope result di- 
rectly from Gaussian fits to the reconstructed Bp distributions such as the one 
shown in Fig. 3. The resulting laboratory momenta are then transformed to 
the projectile system. As it is well known, the centroid values are small (on the 
order of 100-200 MeV/c), i.e. the fragments have practically the same velocity 
as the projectile beam. The central momenta decrease approximately linearly 
with mass loss from the projectile, with a slope similar to what can be found 
in the Morrissey systematics [8]. Since the error bars in our measurements 
are rather large, and do not give more insight than what can be found in the 
literature, we do not discuss this observable here. 



The widths of the longitudinal momentum distributions, c(p\\), are listed in 
Table 2 and are shown graphically in Fig. 5. A marked deviation from the 
parabolic shape predicted by the Morrissey systematics [8] is clearly visible 
above a mass difference of 10 units between fragment and projectile. 
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Table 2 

Longitudinal momentum widths in the projectile system of neutron- deficient frag- 
ments from 36 Ar measured at 1050 A MeV 



Isotope 


<r(p\\) (MeV/c) 


Isotope 


<r(p\\) (MeV/c) 


12 N 


296 ±2 


23 Al 


269 ±6 


13 N 


322 ±3 


23 Mg 


306 ±8 


13Q 


281 ± 10 


24 Al 


293 ±5 


14q 


299 ± 1 


24 Si 


277 ±4 


15q 


289 ±4 


25 Al 


286 ±2 


16Q 


297 ± 16 


25 Si 


270 ±3 


17p 


315 ±4 


26 Si 


253 ±2 


17 Ne 


312 ±3 


27p 


227 ±4 


18p 


288 ± 18 


27 Si 


267 ±3 


18 Ne 


312 ±3 


28 p 


236 ±3 


19 Ne 


286 ±6 


29 p 


216 ±2 


20 Mg 


265 ±4 


29 s 


233 ±2 


20 Na 


323 ±4 


30 S 


205 ±2 


20 Ne 


302 ±8 


31 C1 


165 ±2 


21 Mg 


287 ±8 


31 S 


197 ± 7 


21 Na 


273 ± 16 


32 C1 


182 ±3 


22 Al 


290 ± 1 


33 Ar 


182 ±4 


22 Mg 


305 ±4 


33 C1 


193 ±2 


22 Na 


329 ± 10 


34 Ar 


124 ±3 



4 Discussion 

4-1 Production cross- sections 

The present data constitute an excellent data base to benchmark the predic- 
tive power of different fragmentation models at the limit of unbound nuclei. 
At these energies, the fragmentation process is described in terms of a two- 
step model. In the first step, the collision between the nucleons in the overlap 
region between projectile and target leads to the formation of two residues 
of both reaction partners called prefragments, with similar kinematical prop- 



9 



erties, smaller mass and a given excitation energy. In the second step, the 
residues of both reaction partners release their energy excess by evaporating 
nucleons and 7-rays. The proximity of the investigated nuclei to the proton 
drip line is expected to increase the sensitivity of their production cross sec- 
tions to the excitation energy induced in the first stage of the collision or the 
competition between proton and neutron evaporation during the second stage. 

To investigate the reaction mechanism leading to the formation of the mea- 
sured nuclei we have performed calculations using two different models to de- 
scribe the first stage of the reaction, the intranuclear-cascade model ISABEL 
[13] and the abrasion model of Gaimard and Schmidt [11]. Both models were 
coupled to the ABLA evaporation code [12] to describe the second stage of 
the reaction. These two models follow completely different approaches. In the 
intranuclear-cascade model ISABEL, the interaction between projectile and 
target is described as a sequence of elastic or inelastic nucleon-nucleon colli- 
sions. The moving nucleons follow classical trajectories and the only quantum 
mechanical ingredient is the Pauli blocking. Nucleons with sufficiently high 
momenta can leave their respective nuclear volumes, the leftover nucleons 
form the prefragment. The final prefragment excitation energy is obtained 
from particle-hole excitations in the initial Fermi distribution of the removed 
nucleons plus the energies of the scattered nucleons that end up below the 
"cutoff energy"; the latter has been introduced to terminate the history of 
cascade nucleons with energies too low to escape their respective potential 
wells [13]. 

The abrasion model follows geometrical considerations, the size of the projec- 
tile and target prefragments is given by the non-overlapping region between 
both partners, which is defined by the impact parameter. The excitation en- 
ergy induced in the collision is determined by the energy of the holes of the 
removed nucleons, assuming a Fermi distribution of nucleons inside the nu- 
cleus, which has been empirically demonstrated to be on average around 27 
MeV per abraded nucleon [22] . 

Other than the physical models discussed above, the empirical parametriza- 
tion of fragmentation cross sections, EPAX [7], has little physical input and 
merely aims at reproducing measured data by numerical expressions. Such 
an approach is nevertheless useful e.g. for planning experiments or for cases 
where fast iterative calculations of secondary cross sections are needed. In the 
present context, the EPAX formula is particularly useful since it contains a 
parametrization of the "memory effect" , i.e. the influence of the proton excess 
of the projectile on the proton excess of the fragments. As explained in detail 
in Ref. [7], the numerical form of this "memory effect" has been adjusted to 
measured 58 Ni fragmentation cross sections. A comparison with the present 
high-precision 36 Ar cross sections can serve as a check of the predictive power 
of EPAX for proton-rich projectiles. 
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In figure 4 we compare the measured production cross sections with the re- 
sults of the calculations obtained with the intranuclear-cascade model IS- 
ABEL (dotted lines) and the abrasion model of Gaimard and Schmidt ABRA 
(dashed lines) coupled with the ABLA evaporation model. The predictions 
obtained with the semi-empirical parametrization EPAX are included in this 
figure as well (solid line) . 

In general the three calculations provide a very good description of the mea- 
sured cross sections even for weakly bound nuclei at the proton drip line. Only 
a detailed comparison shows some trends. The ISABEL calculation seems to 
overestimate the production cross sections of heavy-projectile residues while 
the abrasion model provides an excellent description of all the data. Since both 
calculations use the same evaporation model the observed deviations must be 
due to the different models used to describe the first stage of the collision. In 
addition, residues close to the projectile do not reach the limiting fragmenta- 
tion regime and consequently their production cross section is sensitive to the 
entrance channel. The observed overestimation of the production cross section 
with ISABEL can be explained if one considers that this code produces frag- 
ments with 20% less excitation energy on average than the ABRA code. As 
pointed out, the production of heavy projectile residues is very much affected 
by the first stage of the collision. In particular, a lower excitation energy leads 
to shorter evaporation chains and consequently an increase of the production 
cross section of heavy residues. The good description of these data with the 
abrasion-ablation model ABRABLA confirms its predictive power and vali- 
dates the physical assumptions used in this code. 



The overall agreement of the measured production cross sections with the 
EPAX formula [7] is excellent, as can be seen in Fig.4. The present experiment 
reaches the proton drip line up to Z = 13, Al. Even drip-line nuclei are 
produced with cross sections that lie well on the EPAX predictions. This 
confirms observations that were made for another neutron-deficient projectile, 
58 Ni [3]. 

Our new data for 36 Ar (N/Z = 1.0) allow to test the EPAX predictions for 
the production of exotic nuclei over a wide range of project ile-N/Z values 
by comparing them to similar data from 43 Ar (N/Z = 2.39) on 12 C at 222 A 
MeV [23]. We assume that for both incident energies, 222 and 1050 AMeV, the 
cross sections are largely energy-independent, so that EPAX can be applied. 
Fig. 6 demonstrates that at least the proton-rich slopes of the 36 Ar-induced 
distributions and the neutron-rich slopes of the 43 Ar-induced distributions are 
correctly predicted by EPAX. Similarly good agreement was found for isotope 
distributions from 96 Ru and 96 Zr fragmentation [24]. 
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Fig. 6. Production cross sections of P isotopes (left) and S isotopes (right) from the 
fragmentation o/ 36 Ar-f 9 Be at 1050 A MeV (circles, this work) in comparison to 
data from 43 Ar-/- 12 C fragmentation at 222 A MeV (squares, Ref. [23]). The curves 
denote the respective EPAX predictions. 

4-2 Momentum distributions 



A precise knowledge of the kinematics of high-energy projectile fragmentation 
is essential to correctly predict the yields of secondary beams from an in-flight 
separator like the FRS. We have discussed in subsect.3.i how the transmission 
of the FRS is affected by the momentum distributions. 

In the traditional two-step model of high-energy projectile fragmentation, the 
momentum widths are largely determined by the random addition of the Fermi 
momenta of the nucleons that are knocked out in the first ( "abrasion" ) phase 
of the reaction. The corresponding prefragment momentum widths have been 
given by Goldhaber [25] as 



<?{P\\) = °o 



\ 



ApfiAp-A 

An - 



Pf) 



Pf_ 
V5\ 



A pf (A p - A pf ) 
A n -1 



(1) 



where A p (A p f) refers to the projectile (prefragment) mass; the latter cannot 
be observed. The quantity pp denotes the Fermi momentum of the nucleons 
in the projectile. Subsequent evaporation ("ablation") of nucleons in the sec- 
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ond phase additionally broadens the momentum distributions somewhat, but 
the most important effect here is that the number of nucleons is reduced, so 
that the original "Goldhaber" distribution becomes flatter and can be better 
described (at least for fragments close to the projectile mass) by Morrissey's 
formula [8]: 



where now Af denotes the (observed) fragment mass. Systematic studies of 
the validity of the Morrissey systematics have been undertaken about a decade 
ago at the FRS (e.g. Ref. [26]). 

Up to about a mass loss of 10 units, the present data follow nicely the Morrissey 
systematics (dashed curve in Fig. 5). For fragments further removed from the 
projectile, however, the widths tend to saturate. Similar observations were 
made in a comprehensive study of 86 Kr fragmentation at 500 A MeV [4], 
where it was found that the widths even became smaller again for fragments 
with masses smaller than about half the projectile mass. One practical solution 
is to keep the Goldhaber formula (Eq. 1) but fit the width coefficient er to 
the measured data of Fig. 6 (full line). This leads to a parameter of a[ = 
98.2 ± 0.2 MeV/c, smaller than the Goldhaber prediction (for prefragments) 
of = 112 MeV/c. The latter value is based on a Fermi momentum in 
36 Ar of pf = 250 MeV/c [27]. The smaller cq value reflects the evaporation 
mass loss as discussed above. It compares well with other data found in the 
literature for light nuclei, e.g. Viyogi et al. (Ref. [1]) found er = 94 ±5 MeV/c 
for 40 Ar at 213 AMeV. 



5 Summary 

We have studied the fragmentation of neutron-deficient 36 Ar at relativistic 
energies and have obtained precise values for production cross sections and 
longitudinal momentum widths of very neutron-deficient fragments, reaching 
the proton-drip line up to aluminum. Both observables can be reproduced well 
by their respective empirical parametrizations: the cross sections match excel- 
lently those predicted by the EPAX formula, thus confirming its validity also 
for the case of very neutron-deficient projectiles. The longitudinal momentum 
widths are best fitted by a modified Goldhaber formula, where the width con- 
stant er is replaced by an empirical value. The measured production cross 
sections are also compared with physical-model calculations performed with 
(i) the ISABEL intranuclear-cascade model, and (ii) the Gaimard- Schmidt 
abrasion model, both coupled to a standard evaporation code (ABLA). The 




(2) 
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ISABEL code gives an overall good description of the data however, some dis- 
crepancies are observed for the heavier residues. For these nuclide very close 
to the projectile the ISABEL code clearly overestimate the production cross 
sections. This result is interpreted as a consequence of an underestimation of 
the excitation energy induced in the reaction; lower excitation energies lead 
to shorter evaporation chains which increase the production of residues close 
to the projectile. The abrasion model ABRABLA provide an extremely good 
description of all the data. 

These results indicate that in the range of fragment masses studied here the 
assumption of a two-stage model works well; the model consists of a fast 
removal of quasi-free nucleons during the geometrical overlap of projectile 
and target, followed by evaporation of particles from a thermally equilibrated 
prefragment. In addition, we confirm the predictive power of the abrasion code 
that can also be used to describe the production of residual nuclide in more 
complex reactions where the EPAX formula can not be applied, as it is the 
case for fissile projectiles. 
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